Full transmission through perfect-conductor subwavelength hole arrays 
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Light transmission through 2D subwavelength hole arrays in perfect-conductor films is shown to 
be complete (100%) at some resonant wavelengths even for arbitrarily narrow holes. Conversely, 
the reflection on a 2D planar array of non-absorbing scatterers is shown to be complete at some 
wavelengths regardless how weak the scatterers are. These results are proven analytically and 
corroborated by rigorous numerical solution of Maxwell's equations. This work supports the central 
role played by dynamical diffraction during light transmission through subwavelength hole arrays 
and it provides a systematics to analyze more complex geometries and many of the features observed 
in connection with transmission through hole arrays. 

PACS numbers: 42.25.Fx,42.79.Ag,41.20.Jb,78.66.Bz 
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The fraction of light transmitted through a periodic 
array of subwavelength holes perforated in a metallic 
film has been shown to exceed the open fraction oc- 
cupied by the holes at certain wavelengths related to 
the array periodicity and the observed transmis- 
sion can be anomalously large as compared to well es- 
tablished predictions for isolated holes Q. This phe- 
nomenon has triggered an intense amount of activity 
on both experimental El, IS SJSJe , & IS HI an d theo- 
retical 0, i, US EH, El EE 111 H! fronts. While the 



first observations were made in the near infrared (NIR) 
0, 0, H, H, 0| , the effect has been recently corroborated at 
THz frequencies [S IS , where claims have been made 
that the transmission is even larger [S 0] . Two different 
kinds of complementary interpretations of the enhanced 
transmission have been proposed so far, depending on 
whether surface plasmons ElEL 111 ll, 13 , 14 1 or dynam- 
ical diffraction resonances [5|, [lCj, Il5| are invoked as the 
origin of the effect. 

In this Letter, Babinet's principle [r| is used to 
study light transmission through hole arrays in perfect- 
conductor thin films (PCTFs) by relating it to the re- 
flection on planar arrays of metallic disks, which are 
solved from the multipolar polarizability of single disks. 
It is found that transmission reaches 100% at some res- 
onant wavelengths regardless how small the holes are, 
while the complementary system exhibits perfect reflec- 
tion resonances for arbitrarily small disks. These reso- 
nances are associated to divergences in the coherent in- 
teraction among disks/holes, as in the Wood anomalies 
171 ]. A quasi-bound-state (QBS) close to the onset of 
the first propagating diffraction channel (i.e., close to 
a Rayleigh frequency) shows up due to the interaction 
of the scattered field with diffraction modes. The cou- 
pling between the incoming field and the QBS leads to 
transmission/reflection resonances which display charac- 
teristic Fano-line shapes [IS]. This is rigorously proven 



by analytically solving the small-hole limit and also via 
full numerical solution of Maxwell's equations. A typi- 
cal result of our analysis is advanced in Fig. UJ showing 
full transmission at wavelengths A immediately above the 
lattice period a. The theory of hole arrays in PCTFs 
traces back to the works of Eggimann, Collin, and Chen 
[l9| , and McPhedran et al. [2(| . The latter presents com- 
prehensive comparisons with experiments, demonstrating 
the existence of 100% transmission peaks at wavelengths 
immediately above the period. Here, we extend the result 
of full transmission to arbitrarily small holes. 

Babinet's principle flfl ] provides a powerful tool for 
analyzing complex perfect-conductor planar structures. 
For instance, it implies that any self-complementary 2D- 
patterned flat screen (e.g., a chess-board of square aper- 
tures) reflects and transmits exactly half of the unpo- 
larized light coming at normal incidence for any wave- 
length. Also, a simple derivation of Bethe's result for 
a small isolated hole in a metallic screen [2] follows 
from this principle: for incidence normal to a metallic 
disk of diameter D, the electric polarizability reduces to 
olb = D 3 /6ir, while the magnetic response vanishes alto- 
gether [16J, leading to a forward scattering cross section 
given by a/(nD 2 /4) = 4(fcD) 4 / '27 V 2 (normalized to the 
disk area) , as calculated from the Poynting vector in the 
forward hemisphere, where k = 2ir/X; now, Babinet's 
principle implies that this result must coincide with the 
transmission cross section of a hole of the same diameter, 
in perfect agreement with Bethe's formula 0]. 

For hole arrays in PCTFs, Babinet's principle connects 
their transmittance to the reflectance of the complemen- 
tary geometries consisting of planar arrays of perfect- 
conductor thin disks [IS]. In the small- hole limit, this 
can be obtained by self-consistently solving for the elec- 
tric dipoles induced on the disks, pr = a£E lnc (R) + 
Q£ Er'/r^£ R Pr/i where R and R' label disk sites, 
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FIG. 1: Light transmittance through a square array of cir- 
cular holes perforated in a thin perfect-conductor screen as a 
function of wavelength A normalized to the lattice constant 
a for two different diameters of the holes: (a) D — a/2.5 
and (b) D = a/4.5. The light is impinging normal to the 
screen and a 100% transmission peak is obtained in both 
cases. Solid curves: full numerical results. Broken curves: 
analytical model of Eq. ([3]). 



oiE is the electric polarizability of the disks, and G^ R 
is the dipole-dipole interaction dyadic. At normal inci- 
dence on periodic arrays, the polarization is the same in 
all disks, parallel to the incident electric field E mc (|| x), 
and given by 



P : 



l 

OtE 



G f 



where 



Gr = 



E( fc2 

R#0 



+ di x y kR /R 



(i) 



(2) 



Below the threshold of the first diffraction channel, the 
far-field induced by the lattice of identical dipoles is read- 
ily calculated as E md = (2irik/A)p e ife l z l, where A is the 
area of the unit cell and z is chosen perpendicular to 
the lattice plane. From here, one obtains the transmit- 
tance |1 + (2-Kik/A)p\ 2 and reflectance \(2Ttk/A)p\ 2 of the 
disk array, which have to sum 1 because perfect con- 
ductor cannot dissipate energy. This condition leads to 
Im{l/a!£; — Ge} = —2-jrk/A, which generalizes the opti- 
cal theorem for non-absorbing particles, lm{ — \/ole} = 



2/c 3 /3. We can then write the reflectance of the disk 
array, which equals the transmittance of the complemen- 
tary hole array from Babinet's principle, as 

1 
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(3) 



the structure 



Like in the Wood anomalies 
factor Ge diverges when one of the reflected beams 
goes grazing, and in particular near A = a, the pe- 
riod for square lattices. More precisely, Re{G,Ea 3 } = 
4-7r 2 \/2/^A/a - 1 + G(A/a), with A > a, where the first 
term is derived analytically from the divergent terms of 
the sum over reciprocal lattice vectors into which Eq. @ 
can be recast, and C(6) « 35 exp[— 22(9 — 1)] — 118 is a 
smooth fit to the remaining non-divergent terms. 

Obviously, Ge diverges as A — > a + , and thus, given 
an arbitrarily-small hole diameter, there is always one 
wavelength A > a for which the second term in the de- 
nominator of Eq. ([3]) vanishes, implying that the trans- 
mittance (reflectance) becomes 100% at that frequency 
for the array of holes (disks), as shown in Fig. Q] (dashed 
curves). This is graphically illustrated in Fig. [21 where 
Re{G_E} has been represented [Fig. E^a)] and compared 
with Re{l/a£;} (horizontal line). The crossing points 
of these two agree very well with the transmission max- 
ima predicted by the analytical Eq. ^ [dashed curve 
in Fig. HJb)]. Our findings support recent experiments 
performed at THz frequencies, where nearly perfect- 
conductor behavior is expected, and up to 90% absolute 
transmission has been observed @]. 

This divergence of Ge is the same one that produces 
Wood's anomalies and is related to the geometry of the 
system. For disks, this is simply coming from the ac- 
cumulation of in-phase contributions from long-distance 
coherent multiple scattering (MS) [l(J 15 1. In the 
complementary structure (hole arrays), the same sys- 
tematics is at work via Babinet's principle, as described 
above. Moreover, the reflectance (transmittance) of the 
disk (hole) array [Eq. ©] exhibits the same behavior as 
Fano reflection resonances, in agreement with reported 
descri ptio ns of the mechanism of transmission enhance- 
ment 15]. In general, Fano resonances [l8| occur when 
the energy (frequency) of the incoming wave is tuned to 
the energy (frequency) of a QBS. Our QBS is confined 
to the lattice and partially leaking into the continuum 
of vacuum light states, giving rise to some broadening of 
the transmission peaks, as observed in Fig.[TJ 

For perfect-conductor disks or hole arrays, we can then 
consider the transmission enhancement as a consequence 
of the coupling of the incoming field with a QBS of geo- 
metrical origin. Whether this is or is not denoted surface- 
plasmon constitutes in our opinion a matter of taste re- 
garding terminology. 

We have also carried out a rigorous solution of 
Maxwell's equations that is represented in Fig. [1] by solid 
curves. In our method, we first obtain the multipolar 





).4 0.5 0.6 0.7 0.8 0.9 1.0 1.1 1.2 

% I a 

FIG. 2: (a) Wavelength dependence of Ge [Eq. and de- 
termination of transmission maxima and minima: the solid 
(broken) curve represents the real (imaginary) part of Ge for 
a square lattice under normal incidence. Transmission min- 
ima (see vertical dotted lines) result from divergences of Ge, 
while transmission maxima (vertical dashed lines) are derived 
from the condition that Ke{Gs} equals the inverse of the po- 
larizability of the hole, (b) Transmittance under the same 
conditions as in Fig. [Ha) over a wider wavelength range. 



scattering matrix of a disk of the same diameter as our 
holes using the boundary element method [2ll ] . Then, we 
solve the problem of a lattice of disks, which is the exact 
complementary screen geometry of our hole array, using 
the layer-KKR approach [22j]. Finally, the solution for 
the hole array is related to the solution for the disk array 
using Babinet's principle, and in particular, the trans- 
mittance for s (p) polarization in the former equals the 
reflectance for p (s) polarization in the latter. 

Fig. [1] shows results for two different ratios of the di- 
ameter D to the lattice constant a. The analytical model 
of Eq. explains well the shift of the transmission res- 
onance wavelength towards a and its narrowing with de- 
creasing D/a. Actually, Eq. ([3]) works better for smaller 
holes, since it is based upon the dipolar part of their re- 
sponse, and it should be exact in the D/a — > limit, 
although higher perfectness of hole shape, size, and po- 
sition will be necessary to observe narrower resonances. 

Of course, the derivation of Eq. © fails when more 
than one transmitted beam is present: 100% transmission 
can only occur for A > a in our square lattices. This is 
shown in the transmission maxima of Fig. [3fb) for A < 
a, where non-normal beams take part of the light flux, 
and where Eq. ([3]) (broken curve) is shown to reproduce 
qualitatively our rigorous calculations (solid curve). 

A closer look into the near field is given in Fig. [3] for 
on-resonance [(a) and (b)] and off-resonance [(c) and (d)] 



FIG. 3: (a) and (c): squared electric field in log scale for light 
incident normal to the hole array considered in Fig. []Ja), un- 
der resonance and non-resonance conditions, respectively (see 
arrows in [Fig. HI a)] . The plane of representation is parallel 
both to the screen normal and to one of the lattice vectors, 
and it intersects the centers of a row of holes. The metallic 
screen and its holes are represented by vertical broken lines, 
(b) and (d): same for the disk array defined as the comple- 
mentary screen with respect to that of (a) and (c). 



scenarios in a hole array [(a) and (c)] and in its com- 
plementary disk array [(b) and (d)]. The effect of the 
on-resonance hole array [Fig. [Ha)] extends up to a dis- 
tance of the order of the lattice period a, beyond which 
the electric field strength is quite uniform. For the on- 
resonance disk array [Fig. E^b)], there is an evanescent 
field beyond the plane of the disks and total reflection 
establishes an interference pattern on the incoming-light 
side. The off-resonance scenario exhibits interference of 
incoming and partly-reflected light, as well as partial light 
transmission. It should be noticed that the field near the 
hole array is more intense on resonance (this is a sig- 
nature of a QBS) as a result of stronger MS, and this 
anticipates stronger absorption if lossy materials are em- 
ployed rather than perfect conductors, in agreement with 
recent experiments 0j- 

When light is coming under oblique incidence, the 
above model has to be supplemented by adding a phase 
exp{ikj" c • R} to the dipoles induced on the different 
hole sites R. And more importantly, both magnetic 
and electric responses will play a role. In fact, any 
perfect-conductor planar object of finite extension such 
as our disks can sustain induced currents and charges 
only within the plane where it is contained, leading to 
induced electric dipoles parallel to the plane and induced 
magnetic dipoles perpendicular to it in the long wave- 
length limit. The latter can only be excited under oblique 
incidence. 

In particular, the magnetic polarizability of a disk, 
UM = — D 3 /I2ir, reacts only to normal magnetic field 
components. The disks in a lattice are coupled via 
Gm = Er^o exp{ikj^ c • R}(fc 2 + d 2 zz )e lkR /R to produce 
self-consistent perpendicular magnetic dipoles given by 
the magnetic counterpart of Eq. (fTJ). Therefore, there 
are two uncoupled sets of resonances in a disk array that 
respond to external (i) perpendicular magnetic fields and 
(ii) parallel electric fields, which in virtue of Babinet's 
principle are translated into hole-array resonances ex- 
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FIG. 4: Angular dependence of the transmission of s- and 
p-polarized light [region above the light line in (a) and (b). 
respectively] incident upon a perfect-conductor screen perfo- 
rated by an array of circular holes of diameter D = a/2.5 
arranged in a square lattice of constant a. Resonances out- 
side the light cone are also shown by plotting the transmitted 
far field when the screen is exposed to an incident evanescent 
wave. 



cited by external (i) perpendicular electric fields and (ii) 
parallel magnetic fields, respectively. For s-polarized in- 
cident light (E lnc parallel to the screen) only the latter 
can be excited, whereas p-polarized light can couple to 
both types of resonances, as illustrated in Fig. |4] via the 
dependence of the transmittance on incident-light wave- 
length and parallel momentum. The region above the 
light line (A < fc|i nc ) shows transmission maxima (bright 
areas) and minima (right above the maxima), with 100% 
transmission achieved only for p polarization and A > a. 
As expected, p-polarized light couples to more resonances 
[of types (i) and (ii)] than s-polarized light [only reso- 
nances of type (ii)] . This conclusions might be relevant to 
explain recent measurements of polarization- and angle- 
resolved transmission in thick arrays of real metal Q with 
similar qualitative behavior as that of Fig. 0J 

Our methods can be trivially extended to other geome- 
tries like finite arrays [|[ and more complex hole shapes 
0, lH- In particular, the asymmetry of the polarizabil- 
ity tensor of a rectangular metallic plate, which is larger 
along its long axis, is consistent with the stronger red 
shift of the transmitted maxima observed for lattices of 
square holes Q as compared to circular holes of similar 
size. Also, new metamaterials composed of layers of hole 
arrays and disk arrays could be analyzed by straightfor- 



ward generalization of the methods presented above. 

In real metals, the dispersion relation of surface plas- 
mons deviates from the light line and it is plausible that 
the transmission enhancement observed in the NIR is the 
result of in-phase MS mediated by electromagnetic prop- 
agation via these modes. 

Furthermore, the metal films are supported on glass 
substrates in many experiments [3, 0, 0] ■ This analysis 
can be easily extended to a PCTF surrounded by di- 
electrics of different index of refraction on either side as 
shown by Dawes et al. [23], in which case full transmis- 
sion resonances occur only when the wavelength in the 
higher-index region is above the smallest reciprocal sur- 
face lattice vector of the array (e. g., near A ~ 1.4a for a 
square lattice supported on glass) , and the effective polar- 
izability of the hole differs from the symmetric case con- 
sidered above, but apart from this the long-wavelength 
analysis presented above remains the same. 

In summary, we have shown that the extraordinary 
transmission of light through hole arrays in PCTFs can 
be fully explained in simple terms by invoking Babinct's 
principle. The existence of 100% transmission maxima 
has been established, regardless how small the holes are. 
The origin of full transmission can be found in phase 
accumulation during long-range dipole-dipole interaction 
among holes when the wavelength is close to the onset 
of of the first diffraction order. Potential application to 
narrow filters can be envisioned in the microwave and 
THz regimes, since the transmission resonances are in- 
creasingly shaper as the holes become smaller. 
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